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Abstract 

The recent development of a statistical model for incompressible Navier-Stokes (NS) fluids based 
on inverse kinetic theory (IKT, 2004-2008) poses the problem of searching for particular realizations 
of the theory which may be relevant for the statistical description of turbulence and in particular 
for the so-called homogeneous, isotropic and stationary turbulence (HIST). Here the problem is set 
in terms of the 1— point velocity probability density function (PDF) which determines a complete 
IKT-statistical model for NS fluids. This raises the interesting question of identifying the statistical 
assumptions under which a Gaussian PDF can be achieved in such a context. In this paper it is 
proven that for the IKT statistical model, HIST requires necessarily that /i must be SIED (namely 
stationary^ isotropic and everywhere- defined). This implies, in turn, that the functional form of the 
PDF is uniquely prescribed at all times. In particular, it is found that necessarily the PDF must 
coincide with an isotropic Gaussian distribution. The conclusion is relevant for the investigation 
of the so-called homogenous, isotropic and stationary turbulence. 

PACS numbers: 05.20Jj,05.20.Dd,05.70.-a 



I. INTRODUCTION: HYDRODYNAMIC TURBULENCE 



The turbulence problem of hydrodynamics is made difficult by the fact that there does 
not exist a definition of the turbulent state that is universally accepted. The concept of 
fluid state is nevertheless well defined in fluid dynamics, as it follows from the theory of con- 
tinuum media. In particular, it is prescribed by an appropriate set of suitably smooth real 
functions {Z} = {Zi,i = l,n} denoted as fluid fields, which must be physically realizable, 
i.e., identified with physical observables. The fiuid fields by assumption satisfy a well-posed 
initial-boundary value problem represented by a set of PDEs denoted as fiuid equations. 
In the theory of continua this means that they are necessarily described by strong solu- 
tions (i.e., they are defined and are least continuous everywhere in the existence domain). 
The same type of requirement is also manifestly imposed by comparison with experimental 
observations. 

Despite these premises, a precise mathematical definition of the concept of "turbulent 
state" is still missing. In experiments on isolated incompressible viscous fiuids, turbulence is 
typically associated with the manifestation of fiuid motion in which the state of the fiuid is 
"turbulent" , i.e., is in some sense random and decays in time until the fiuid comes ultimately 
to a state of rest. In the past at least three different views have been adopted in this regard. 

The first one is the mathematical viewpoint represented by the Leray theory of turbulence 
(see Stewart, 1988 l|) according to which "...turbulence is a fundamentally different problem 
from smooth flow" , which cannot be described in term of strong solutions. This conjecture 
actually lead Leray (1931, |2"4|) to introduce the concept of weak solutions, i.e., solutions 
which are not deflned, and are not continuous, everywhere in the existence domain. Although 
appealing and extremely fruitful for its mathematical implications the theory is manifestly 
un-physical and therefore should be rejected. 

The second one is the so-called deterministic theory of turbulence, for which turbulence 
should be produced, instead, by the occurrence of a strange attractor (Ruelle-Takens theory; 
Berge, 1984 jsj]). According to this view, there should exist a classical dynamical system, 
characterizing in some sense the time evolution of the fluid flelds, which should exhibit a 
chaotic behavior on a suitable invariant (hyper-)surface. The basic manifestation of chaos 
for such a dynamical system would be the occurrence of an inflnite cascade of quasi-periodic 
phenomena in which inflnitely many periods are sequentially generated in order to give the 
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appearance of randomness (Hopf-Landau theory; Landau 1959 [6|). 

However, the third, and most popular approach, is probably the statistical theory of 
turbulence, which historically can be referred primarily to the work of Kolmogorov (Kol- 
mogorov, 1941 and Hopf (Hopf, 1950/51 10|). The statistical treatment of fluids usually 
adopted for turbulent flows (which may be invoked, however, to describe also regular flows) 
consists, instead, in the introduction of appropriate axiomatic approaches denoted as statis- 
tical models {/, F}. 

Their construction involves, besides the speciflcation of the phase space (F) and the 
probability density function (PDF) /, the identiflcation of the functional class to which / 
must belong, denoted as {/} . As a consequence, the complete set of fluid flelds {Z} , or 



only a proper subset as in the case of the Monin-Lundgren approach 
in terms of suitable functionals (called moments) of /. 
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are expressed 



A. Stochastic representation of turbulence 



A basic aspect of the statistical description is the introduction of an ensemble average 
operator acting on the possible realizations of the fluid. A possible deflnition of such operator 
may be achieved by representing the fluid flelds in terms of hidden variables [l^, lisl. By 
deflnition they denote a suitable set of independent variables a = = 1, fc} G C R*^, 
with k > 1, which cannot be known deterministically, i.e., are not observable. In the context 
of turbulence theory these variables are necessarily stochastic. This means that they are 
characterized by a suitable stochastic probability density g deflned on (see definitions and 
related discussion in the Appendix, Subsection B), while the ensemble average (■) can be 
identified with the stochastic-averaging (■)^ defined by Eq. flSOj) [see Appendix A]. Hence, 
for turbulent fiows the fiuid fields - together with the PDF/iand the vector field F(x, t, a) 



can be assumed to admit a representation of the form 
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19| 



{Z} = {Z{r,t,a)} 

/i = /i(r, u,t,a) (1) 
F = F(x,t,a) 

to be defined in terms of a set of hidden variables a and a stochastic model {g,Va} (see 
again Subsection B in the Appendix). Hence, {Z} , fi and F(x, t, a) are necessarily non- 
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observable. Nevertheless, if we assume that the fluid flelds {Z} are uniquely-prescribed 
ordinary functions of (x,t, a) deflned for all (x,t, a) G F x / x Va, it follows that they can 
still be considered conditional observables (see Appendix, Subsection A). Similar conclusions 
apply to /i, and to the vector fleld F(x, t, a) as well. 



B. The Navier-Stokes dynamical system 



It is, however, generally agreed that a common important property should characterize 
all NS fluids, either regular or turbulent: this is related to the existence of the phase- 
space dynamical system, to be denoted as NS dynamical system, which advances in time the 
complete set of fluid flelds deflning the fluid state. 

In other words, should such a dynamical system actually exist, it would permit to cast the 
complete set of fluid equations in terms of an equivalent (and possibly inflnite) set of ordinary 
differential equations which deflne the dynamical system itself. For contemporary science the 
determination of such a dynamical system represents not simply an intellectual challenge, 
but a fundamental prerequisite for the proper formulation of all phenomenological theories 
which are based on the description of these fluids, and hence both for the deterministic and 
statistical approaches to turbulence. These involve^ 



phase-spac_e 
dynamics 
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br example, the understanding of the 
agrangian description of fluids 13|, [ij] relevant to determine tracer-particle 
16| as well as the time-evolution of scalar and tensor flelds in turbulent flows. 



the search of exact (or approximate) kinetic closure conditions for statistical models (such as 
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12|). the investigation of stochastic models 



in the case of the Monin-Lundgren hierarchy 

17,[l8| able to reproduce phenomenological data (such as the two-point velocity increments 
PDFs , the theoretical prediction of multi-point velocity probability densities, all 

essential ingredients in fluid dynamics and in applied sciences. 

Surprisingly, although phase-space descriptions of incompressible fluids described by the 
incompressible NS equations (INSE) have been around for a long time, starting from the 
historical work of Hopf (see also Hop f, 1952 [l^), Edwards (Edwards, 1964 j26|) and Rosen 
(Rosen, 1971 2^), until recently 24 1 the problem [of the search of the NS dynamical system] 

las remained unsolved. Its solution for the incompressible NS equations (see also Refs. 
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25| for its extension to quantum and magneto fluids) is based on the construction of 



a statistical model {/i,F} for the 1-point PDF fi which is required to obey a Liouville 
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equation and whose moments determine - via suitable velocity-moments - 
of fluid fields which describe the state of the fluid. As indicated elsewhere 
can be extended also to the statistical treatment of turbulence theory. 



the complete set 
17] , the approach 



II. MOTIVATIONS 

An unsolved problem in the statistical theory of turbulence concerns the so-called homo- 
geneous, isotropic and stationary turbulence (HIST) arising in incompressible Navier-Stokes 
(NS) fluids. This concerns in particular the determination of the form of the 1— point proba- 
bility density function (PDF), /i, occurring in the presence of HIST, which should uniquely 
determine, in turn, the statistical model {/, T}. 



According to some authors (see in particular Batchelor 20]) this is predicted as almost- 



Gaussian, while others 
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221] have pointed that the tails of the PDF might exhibit a 
strongly non- Gaussian behavior. Despite recent attempts at possible a theoretical explana- 
tion [23], still missing is a definite answer to the question whether a generalized behavior 
of this type should actually be expected or not. Apart insufficient experimental evidence, a 
major difficulty is represented by the lack of a consistent theoretical description of the PDF 
in the presence of HIST, permitting a rigorous definite answer to this question. This raises 
the interesting question whether, in some suitable setting, i.e., for appropriate statistical 
models, the problem can actually be solved. Being a subject of major importance - not only 
in fluid dynamics but also in statistical mechanics and, as we intend to prove, in kinetic 
theory - the issue deserves a careful investigation. The goal of this paper is to pose the 
problem in the framework of the complete inverse kinetic theory (IKT) approach developed 
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19|). 



by Tessarotto et al. [24] for incompressible NS fluids (see in particular 

Here we intend to prove that based solely on IKT the appropriate form of the 1-point 
PDF in the presence of HIST can actually be uniquely established for these fluids, based on 
suitable statistical assumptions stemming from the requirement of existence of HIST and 
appropriate initial conditions and the requirement that the initial 1-point PDF is deter- 



mined imposing PEM (principle of entropy maximization, Jaynes, 1957 SOj; see also related 



discussion in Ref.[25l]). This is shown to be described by a probability density /i, deflned on 
the restricted phase-space T = Q x U, with Q and U denoting respectively the conflguration 
space of the fluid [to be identifled with a bounded subset of the Euclidean space R'^] and the 
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Euclidean velocity space U C R^, which fulfills the following properties (#1-8): 

1. (Property #1) it depends explicitly on the fluid fields {Z{v,t,a)} = 
{V{v,t,a),pi{v,t,a)} , namely is of the type /i = /i(r, v,i, a), with V(r, i,Q;) de- 
noting the fiuid velocity and pi{v,t,a) the kinetic pressure. In particular pi(r, a) 
is defined as the strictly positive function pi(r, a) — p{r, t, a) + po{t, a) + 0(r, t, a), 
with p{r,t,a),po(t,a) and 0(r, t,Q;) representing respectively the fiuid pressure, the 
(strictly-positive) pseudo-pressure and the (possible) potential associated to the con- 
servative volume force density acting on the fiuid; 

2. (Property #2) /i is a velocity probability density, i.e., it is normalized in velocity space 
so that 

j d\Mr,v,t,a) = l; (2) 

u 

3. (Property #3) /i is strictly positive in the velocity space U ; 

4. (Property ^4) fi Galilei invariant in velocity space, namely it is invariant with respect 
to a transformation of the form: 

" (3) 

with Vo such that v -|- Vq e 6'^. As a consequence, fi is of the form 

/i = /i(r,u,t,a), (4) 

with u = u(r, t, a)—v—Y{r, t, a) denoting the relative velocity, namely is homogeneous 
in the velocity space U; 

5. (Property #5) the velocity space U coincides with R^. fi fulfilling this property is said 
everywhere defined (in U = 



6. (Property #6) /i is stationary, namely it can depend on time only via the fiuid fields: 

/i = /i(r,u,a), (5) 

7. (Property #7) /i is isotropic in velocity space, i.e., it is of the form 

/i = /i(r, |u| ,«); (6) 
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8. (Property 7^8) /i is a Gaussian distribution of the form 

/A/(r, |u| ,pi{r,t,a)) = ^ — -exp J — — - )■ , (7) 

a) [ v^hpir,t,a) 

with Vthpijiti ot) = ^y2pl{r, t, a) / po denoting the thermal velocity associated to the 
kinetic pressure. 



III. HOMOGENEOUS, ISOTROPIC AND STATIONARY TURBULENCE 

In fluid dynamics two types of descriptions of the fluid, respectively denoted deterministic 
and stochastic, can be distinguished, in which the fluid fields describing the state of the fluid 
are treated respectively as deterministic or stochastic functions. In both cases the fluid 
fields {Z} = {Zi,i = l,n} are considered, in a suitable existence domain, suitably smooth 
strong solutions of the fluid equations. In the so-called statistical theory of turbulence, 
historically referred primarily to the work of Kolmogorov (Kolmogorov, 1941 (tI) and Hopf 
(Hopf, 1950/51 [lO'i]), turbulence is intended as the characteristic property of the fluid in 
which the fluid fields {Z} can only be prescribed in a statistical sense. This implies that 
they are necessarily stochastic functions of the form {Z} = {Z{r,t,a)} characterized by a 



stochastic PDF g{r,t, a) defined for all {r,t) G fix J [17|,ll9| and a = {oj, i = 1, A;} G K« C R' 



(with k > 1), with a suitable stochastic variables independent of r,t. Thus, introducing 
the stochastic- averaging operator (■)^ = J dag{T,t,a)-, acting on an arbitrary integrable 

function it follows that the fluid fields Zi{r,t,a) (for i = l,n) can always be represented 
in terms of their stochastic decompositions Zi = (Zi)^ + 6Zi, with (Zi)^ denoting their 
stochastic averages. 

A widespread conjecture is that turbulence, at least in special circumstances [usually 
ascribed to the so-called "fully developed" turbulence (FDT)], should be characterized by 

ar, the concept of homogeneous. 



a 



(see also Refs. 
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29|) is related 



certain universal properties. These concern, in particu 
isotropic and stationary turbulence (HIST). Its definition 
to the assumed properties of the operator (■.)^ and of the velocity increments 

dVi = V,{n,t)-Vi{r,t), (8) 

which are assumed to be defined for arbitrary displacements 

dr = ri - r (9) 



such that both r and ri belong to 
1. Definition - HIST 

Turbulence is said homogeneous, isotropic and stationary if: 

• HIST Requirement #1: the stochastic-averaging operator (Zi)^ commutes with all the 
differential operators appearing in the fluid equations (namely, for the NS equations, 
this means it must commute with the operators ^, V and V^); 

• HIST Requirement #2: for all n E Nq the structure functions - i.e., the stochastic- 
averages of {dVi}"' , with sl"'\r,dr,t) = {{dVi}"')^ - are respectively: 

2a) independent of r, namely for alH = 1, 2, 3 

^f) = 5f)(dr,t) (10) 

{homogeneous turbulence) ; 

2f,) independent of the directions of dr and V, hence for alH = 1, 2, 3 : 

5j") = 5W(r,/,t) (11) 

[isotropic turbulence) , where / = |dr| is the magnitude of the displacement ([9]); 
2c) independent of t, hence for alH = 1, 2, 3 : 

^^) = ^f)(r,dr) (12) 

[stationary turbulence). 

Thus, for all n G Nq and i = 1,2, 3, HIST is by assumption characterized by structure 
functions of the form 

5'f)=^W(/), (13) 
i.e., depending solely on the magnitude of the displacement ([9]). 
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IV. THE IKT-STATISTICAL MODEL FOR TURBULENT FLUIDS 



Starting point for the statistical treatment of turbulence in NS fluids in the IKT approach 
is the introduction of a statistical model {/i, F} {to be denoted as IKT- statistical model) for 
the INSE problem. The corresponding fluid fields are {po,y^{r,t,a),p{r,t,a), St} with 
Po and St to be identified, respectively, with the constant mass density and the constant 
thermodynamic entropy. In the following the fluid fields are required to be: 1) they are 
strong solutions of the INSE problem in Q x I x Va, with bounded configuration space Q 
{internal domain) ; 2) global solutions, i.e., defined for alH G / = M. 



The construction of {/i,r} 17|] involves the definition of a suitable PDF /i defined 
on a phase space F, denoted as 1-point velocity PDF, which permits the representation, 
via a suitable mapping {/i,F} =^ {Z} , {Z} denoting the complete set of the fluid fields 
{Z} = {Zi, i = l,n} defining the state of the fluid. In particular by assumption F is identified 
with the restricted phase-space T = Q x U xVa [with closure T = Q x U x V^]; furthermore, 
the 1-point velocity PDF: 

1. is taken of the general form fi{t) = /i(x, t,Q;), with x = (r, v), where r ^ Q and 
V G f/ C R3 (with U defined as the open subset of spanned by v on which /i > 0). 
In addition /i is by assumption Galilei invariant and hence invariant w.r. to (jS]). It 
follows that /i is necessarily homogeneous in velocity space, namely of the form (j3]); 

2. determines in terms of suitable moments the complete set of the fluid fields {Z} which 
define the state of the fluid. This requires that the fluid fields are determined by the 
velocity moments 

j d^vG{x,t,a)fi{x,t,a) =Y{T,t,a),pi{r,t,a), (14) 
u 

defined respectively for the weight-functions G{'x.,t,a) = v,poM^/3, whereas St is 
identified with the Boltzmann-Shannon entropy, i.e., the phase-space moment 

S{h{t)) = - j rfxA(x,t,a)ln/i(x,t,a), (15) 
r 



by imposing that Vt G / = M the constraint 25] 



St = S{h{t)). (16) 
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The set of equations (HM and (fT5|) . denoted as correspondence principle, are assumed do 
be fulfilled identically in the whole existence domain of the fluid fields V(r, t, a),pi(r, t, a) 
and St- This implies manifestly that /i(x, t, a) must be defined and strictly positive on 
r = Q X U xVa ■ The time-evolution of /i(x, is then uniquely determined by the flow 
{stochastic N-S dynamical system) 

Tt,,t : Xo x(t) = Tt„,tXo, (17) 

with Tt^^t the corresponding evolution operator generated by the initial value problem 

^x = X(x,t,a) (18) 

In particular, denoting by J{t) = = exp < J dt'^^jTr ■ X(x(t'), t', a) > the Jacobian de- 



to 



terminant of the flow the time-advanced PDF /i(x(t), t, a) is uniquely determined requiring 
that for all (xo, <y,to, t, ) G F x / x / it satisfies the Lagrangian inverse kinetic equation (IKE) 

/i(t) = J(t)/i(to). (19) 



In particular, as shown elsewhere [19|, the initial PDF fi{to) = f{x.,to,a) can be assumed 
to be a strictly positive smooth PDF of the form 

with h{to) = h{yio,to,a) to be determined and (/i(to))n subject to the constraint 

{Mto))n = ff'''\t^.^o,a) (21) 

{physical realizability condition). Here /[■^^^'^\to,Vo,C() denotes a suitable continuous 
velocity-frequency density, uniquely associated to the initial fluid velocity Y{ro,to,a). Fur- 
thermore, {fi{t))Q denotes the configuration-space average (defined at at time t), i.e., 

(/i(r,v,t,a))^ = j dhf,{T,v,t,a), (22) 
while fi{Q) = J d^r > is the canonical measure of Q. As a consequence h{to), can be 



uniquely determined imposing that /i(xo, to, «) satisfies PEM 30|], namely the variational 
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principle SS{fi(t)) = subject to the constrains (l20l) and (12T|) . taking the form of a gener- 
alized Gaussian distribution 



h{to) = exp [-1 - Xo{ro,to,a) - X2i{ro,to,<y)u^{ro,to,<y)} ■ (23) 

Here, \o{ro,to, a) and \2i{ro,to, a) (for i = 1,2, 3) denote suitable Lagrange multipliers to be 
determined imposing the moment equations f|T^ and f|T5l) . together with the constraint f|T6l) . 
Furthermore, /i(to) can be shown to take the form of an isotropic Gaussian distribution, i.e., 

ho{t) = exp { — 1 — Ao(r,t, a) — A2(r,t, a)'u^(r,t, a) } (24) 

[see related discussion the Appendix]. 

The equivalence theorem pointed out elsewhere 3l| between the INSE problem and the 
NS dynamical system ( |T7|) then warrants the validity of the correspondence principle, i.e., 
that the moments prescribed by equations f|T^ and f|T5l) actually define a strong solution of 
the INSE problem. 



V. CONNECTION WITH HIST - SIED IKT-STATISTICAL MODELS 

Here we are interested in determining a particular subclass of IKT-statistical models 
{/i,r} which may fulfill at least some of the properties which characterize HIST for NS 
fluids. Let us now analyze the consequences placed by (the assumption) of HIST. First, let 
us require that - consistent with the requirement of FDT - the velocity space U [on which 
/i is defined and strictly positive] coincides with and hence /i is everywhere defined. 
Second, we notice that the constraints imposed by Eq.(lT3l) on arbitrary structure functions 
S'j-"'' can generally be satisfied only if: 

• flit) if stationary in the sense ([5]); 

• /i (t) isotropic in velocity space (since no preferred direction in velocity space can exist 
in such a case). 

As a consequence, in the presence of HIST, /i is necessarily stationary, isotropic and 
everywhere defined (SIED). In the following we shall denote as SIED the IKT-statistical 
models {/i,r} which fulfill these requirements. In such a case it is immediate to reach the 
following result [which proves Properties #1-8]: 
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THM.l - Characteristic property of SIED {/i,r} 

For SIED IKT- statistical models {/i,r} there results identically in T x I : 



Kit) = /Af(r, |u(r,t,a)| ,pi(r,t,a)), (25) 
1, (26) 



where /M(r, |u(r,t,a)| ,pi{r,t,a)) is the Gaussian 1-point PDF ^ and ho{t) the Gaussian 
distribution [24\ )- In particular, X{r,t,a) and X2{r,t,a) (for i = 1,2,3) denote suitable 
Lagrange multipliers to be determined imposing the moment equations [l4\ ) and l[TS\) . together 
with the constraint [T^) . 

PROOF In fact, first, since fi{t) is stationary for all t G / = M, fi{t) is necessarily deter- 
mined by PEM , subject to the constrains (120|) and fl2T]) . Therefore for all t G / = R, /i(t) 
has the form 

with h{t) = h{v,w,t,a) of the type (!23l) . Furthermore, since {/i,r} is by assumption 
isotropic and everywhere-defined (so that U = M^) it follows necessarily for all (v,t) G 
]R3 X / that 



c{t,a)>0, (28) 



{Ht))n 

h{t) = exp {-1 - Ao(r,t, a) - A2(r,t, a)u^} = ho{t). (29) 

Here c{t, a), Xoij.t, a) and A2(r,t, a) denote, respectively, a suitable strictly positive function 
of time and two Lagrange multipliers defined so that there results identically in x / 



d'vh{t) = 1, (30) 
j d\P^h{t)=p,{T,t,a), (31) 

R3 

while the pseudo-pressure po{t,a) must be determined by imposing the constraint f|T6l) . 
Hence, it is always possible to set c{t,a) = 1, so that (!26l) is identically satisfied. As a 
consequence, it follows identically that ho{t) = /^(r, |u(r,t, a)| ,pi(r, t, a)), which proves 
Eg. (123!) too. Q.E.D. 
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VI. CONCLUSIONS 



In this paper properties of the IKT-statistical model, {/i,r}, defined in terms of the 
1-point PDF fi for a turbulent fluid obeying the INSE problem, have been investigated. 
In particular we have proven that for the IKT statistical model {/i, F} : 

• the assumption of HIST requires that /i must necessarily be SIED, i.e., stationary and 
isotropic and everywhere-deflned in velocity space; 

• the requirement of stationarity imphes that PEM must hold identically for all t e 
7 = R. As a consequence the functional form of the function h{t) remains uniquely 
determined. 

Main result is the proof here achieved (THM.l) that the requirement of fi to be SIED 

implies necessarily that fi must coincides identically with an isotropic Gaussian distribution. 

The conclusion is relevant for the investigation of the so-called homogenous, isotropic and 
stationary turbulence. 
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VIII. APPENDIX: THE MATHEMATICAL DESCRIPTION OF INCOMPRESS- 
IBLE NS FLUIDS 

In fluid dynamics the state of an arbitrary fluid system is assumed to be deflned every- 
where in a suitable extended configuration domain Q x 7 [Q denoting the conflguration 
space and 7 C M the time axis] by an appropriate set of suitably smooth functions {^} , 
denoted as fluid fields, and by a well-posed set of PDEs, denoted as fluid equations, of which 
the former are solutions. The fluid fields are by assumption functions of the observables 
(r,i), with r and t spanning respectively the sets Q and 7, namely smooth real functions. 
Therefore, they are also strong solutions of the fiuid equations. In particular, this means 
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that they are are required to be at least continuous in all points of the closed set Qx I, with 
= ^7 U dQ closure of Q. In the remainder we shall require, for definiteness, that: 

1. Q [configuration domain) is a bounded subset of the Euclidean space on K , 

2. I [time axis) is identified, when appropriate, either with a bounded interval, i.e., 
/=]to,^i[ ^ IR, or with the real axis M; 

3. in the open set Qxl the functions {Z} , are assumed to be solutions of a closed set 
of fluid equations. In the case of an incompressible Navier-Stokes fluid the fluid fields 
are {Z} = {V,p, St} and their fluid equations 

P = Po, (32) 

V ■ V = 0, (33) 

iVV = 0, (34) 

= 0, (35) 

Z{r,to) = Z,{r), (36) 

^(i",i)lan = ^^(r,i)lan' (37) 

where Eqs. (!32|) - (!35|) denote the incompressible Navier-Stokes equations (INSE) and 
Eqs. fl32|) - fl37|) the corresponding initial- boundary value INSE problem. In particular, 
Eqs. (132|) - (!37j) are respectively the incompressibility, isochoricity, Navier-Stokes and 
constant thermodynamic entropy equations and the initial and Dirichlet boundary 
conditions for {Z} , with {Zoir)} and {Z^„(r,t)|^j^} suitably prescribed initial and 
boundary-value fluid fields, defined respectively at the initial time t = to and on the 
boundary dfl. 

4. by assumption, these equations together with appropriate initial and boundary condi- 
tions are required to define a well-posed problem with unique strong solution defined 
everywhere in Qxl. 

Here the notation as follows. is the NS nonlinear operator 

Ny = ^V-FH, (38) 
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with V and Fh denoting respectively the Lagrangian fluid acceleration and the total force 
per unit mass 

^V=|v + V.VV. (39) 

¥h = - — Vp + — f + vV^Y, (40) 

Po Po 

while Po > and z/ > are the constant mass density and the constant kinematic viscosity. 
In particular, f is the volume force density acting on the fluid, namely which is assumed of 
the form 

f =-V0(r,t) + f«, (41) 

4>{r,t) being a suitable scalar potential, so that the first two force terms [in Eq. (HO|) ] can be 
represented as — Vp + f = — Vpr + ^r, with 

p,(r,t) =p(r,t)-0(r,t), (42) 

denoting the reduced fluid pressure. As a consequence of Eqs.(l32l).(l33l) and ( !34|) it follows 
that the fluid pressure necessarily satisfies the Poisson equation 

V'p = S, (43) 

where the source term S reads 

S = -p„V ■ (V ■ VV) + V ■ f . (44) 
A. Physical/conditional observables - Hidden variables 

The fluid fields {Z} are, by assumption, prescribed smooth real functions of (r, t) e 
Q X I. In particular, they can be either physical observables or conditional observable, 
according to the definitions indicated below. 

1. Definition - Physical observable/ conditional observable 

A physical observable is an arbitrary real-valued and uniquely-defined smooth real func- 
tion of (r,t) G r2 X J. Hence, as a particular case (r,t) are observable too. 
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A conditional observable is, instead, an arbitrary real-valued and uniquely-defined smooth 
real function of (r,t) E I which depends also on non-observable variables and is, as such, 
an uniquely-prescribed function of the latter ones. 



Therefore the functions Zi can be assumed respectively of the form 
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(45) 



or 



(46) 



« G ^ IR'^ (with A; > 1) denoting a suitable set of hidden variables. In fluid dynamics 
these are intended as: 



2. Definition - Hidden variables 

A hidden variable is as an arbitrary real variable which is independent of (r, t) and is not 
an observable. 

B. Deterministic and stochastic fluid flelds 

Hence, fluid fields of the type ( H6|) are manifestly non-observables. However, if in the 
whole set Q X I X V^, they are uniquely-prescribed functions of (r, t,a) then they are con- 
ditional observables. Hidden variables can be considered in principle either deterministic or 
as stochastic variables, in the sense specified as follows. 

1. Definition - Stochastic variables 

Let {S, S, P) be a probability space; a measurable function a iS — > Va, where ^ 1^^, 
is called stochastic (or random) variable. 

A stochastic variable a is called continuous if it is endowed with a stochastic model 
{da, Va} , namely a real function ga (called as stochastic PDF) defined on the set Va and 
such that: 

1) Qa is measurable, non-negative, and of the form 

9a = 9a{r,t,-); (47) 
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2) if A C is an arbitrary Borelian subset of Va (written A 6 BiVa)), the integral 



Pa{A) = j rfx5(„(r,t,x) 



(48) 



exists and is the probabihty that a Gv4; in particular, since a G V^, (^q, admits the normal- 
ization 

(49) 



j dxgair,t,x) = Pa(\4) = 1. 



The set function : BiVa) — )■ [0, 1] defined by fH5]) is a probabihty measure and is called 
distribution (or law) of a. Consequently, if a function f:Va — > Vf C is measurable, / 
is a stochastic variable too. 



Finally define the stochastic- averaging operator (■)^(see also 17|, ll9|) as 



(/L = (/(y,-)L = J rfx(7„(r,t,x)/(y,x) 



(50) 



for any Pa-integrable function /(y, ■) : — M, where the vector y is some parameter. 



2. Definition - Homogeneous, stationary stochastic model 

The stochastic model {ga, Va} is denoted: 

a) homogeneous if ga is independent of r, namely 

ga = gait,-); (51) 

b) stationary if ga is independent of t, i.e., 

ga = 9a{rr)- (52) 



3. Definition - Deterministic variables 

Instead, if ga{r,t, ■) is a deterministic PDF, namely it is of the form 

(7„(r,t,x) = (5(^)(x-a,), (53) 

^{*:)(^x — Oo) denoting the A;- dimensional Dirac delta in the space Va, the hidden variables 
a are denoted as deterministic. 
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Let us now assume that, for a suitable stochastic model {ga,Va}, with ga non- 
deterministic, the stochastic variables Zi = Zi{r,t,a) and /i(r,v,t, a) (where Zj(r, t, ■) 
and /i(r, v,t, ■) are measurable functions) admit everywhere in Q x / and T x I the 
stochastic averages {Zi) ^ and (/i)^ defined by ( l50l) . 

Hence, Zi = Zi{r,t,a), /i(r, v,t, a) and the mean-field force F(/i) [see Sections 2,3 and 
4] admit also the stochastic decompositions 

Z, = {Z,)^ + 6Z,, (54) 
/i = + (55) 
F(/i) = (F(A)), + <5F(/i). (56) 

In particular, unless ga{r,t, ■) is suitably smooth, it follows that generally (Zi)^ ,6Zi and 
respectively (/i)^ , Sfi may belong to different functional classes with respect to the variables 
(r,t). 

C. Deterministic and stochastic INSE problems 

Therefore, assuming, for definiteness, that all the fluid fields Z, the volume force f and 
the initial and boundary conditions, are either deterministic or stochastic variables and both 
belong to the same functional class, i.e., are suitably smooth w.r. to (r, t) and a, Eqs. ( |32l) - 
1371) define respectively a deterministic or stochastic initial-boundary value INSE problem. In 
both cases we shall assume that it admits a strong solution in Q x I [ot Q x I x Va)- 

In the first case, which characterizes flows to be denoted as regular, the fluid fields are 
by assumption physical observables, i.e., uniquely-defined, smooth, real functions of (r,t) e 
Q X I [with Q, the configuration space, and Q its closure, to be assumed subsets of the 
Euclidean space on and /, the time axis, denoting a subset of M]. 

In the second case, characterizing instead turbulent flows, the fluid fields are only condi- 
tional observables (see again Subsection A). In this case, besides (r, t), they may be assumed 
to depend also on a suitable stochastic variable a, (with a gV^ and Va subset of M'^ with 
A; > 1). Hence they are stochastic variables too. 
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